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ABSTRACT 

We present the first study of rapidity gaps in e + e~ annihilations, using Z° decays 
collected by the SLD experiment at SLAC. Our measured rapidity gap spectra fall 
exponentially with increasing gap size over five decades, and we observe no anomalous 
class of events containing large gaps. This supports the interpretation of the large-gap 
events measured in pp and ep collisions in terms of exchange of color-singlet objects. 
The presence of heavy flavors or additional jets does not affect these conclusions. 
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Since the initial observation of hadronic jets rapidity has been used to characterize 
the momentum of particles in jets in a frame-invariant manner. The rapidity distri- 
bution has been studied in e + e~ annihilation, ep and hadron-hadron collisions, and 
fixed-target experiments, and is a characteristic of strong interactions that is well de- 
scribed by perturbative QCD combined with iterative models of jet fragmentation 
Hard scattering of quarks or gluons can be modeled by color fields between the outgo- 
ing partons that fragment into the observed final-state particles, typically populating 
the whole rapidity range. 

Recently, exchange of color-singlet objects in hard diffractive hadron-hadron scat- 
tering processes characterized by events containing large gaps in the particle rapidity 
spectrum has been discussed g. Subsequent studies of pp collisions at the Fermi- 
lab Tevatron collider found that roughly 1% of events comprising at least two high- 
transverse-energy (Et) jets contain a large rapidity region between the two highest-i?T 
jets with no particle activity ||, f|, |5|. These events have been interpreted in terms 
of color-singlet exchange between the interacting partons. Exchange of electroweak 
bosons is estimated to contribute only a small fraction of the observed rate of gap 
events and a model incorporating pomeron exchange is in agreement with the data 
J|, |||. Large rapidity gaps have also been observed in roughly 10% of all photopro- 
duced dijet events in deep- inelastic scattering at the HERA ep collider || [7| and have 
also been interpreted in terms of color-singlet exchange. Models involving either vec- 
tor meson dominance of the exchanged virtual photon or hard diffractive scattering via 
pomerons describe the data || [7| . 

Color exchange processes account successfully for the properties of the majority of 
dijet events and may give rise to large rapidity gaps due to random fluctuations. In both 
the ep and pp cases the interpretation of rapidity-gap events in terms of color-singlet 
exchange depends upon an understanding of this color-exchange background. In the 
pp experiments this was estimated by extrapolation of fits to the particle multiplicity 
distribution in rapidity intervals into the zero-particle, or rapidity gap, region j|, ||. In 
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the ep experiments the background was estimated using both a Monte Carlo simulation 
and an ad hoc parametrisation of color-exchange . In both cases direct measurement 
of the spectrum of rapidity gaps arising in color-exchange jet fragmentation is preferable 
and would clarify the interpretation of the large-gap events. 

Electron-positron annihilation into hadronic final states provides an ideal laboratory 
for study of this issue as, in QCD, it proceeds via creation of a primary quark and 
antiquark connected by a color field which fragments into the observed hadrons. The 
inclusive particle rapidity distribution has a broad plateau centred at zero || . Inclusive 
studies of the local rapidity-density of particles have been performed , with the aim 
of investigating scale-invariant cascade mechanisms in multiparticle production, but no 
previous study has been made of the size of rapidity gaps between adjacent particles. 
Large gaps are expected to occur in e + e~ — > Z° — > qqqq or qqgg events when two 
color-singlet pairs of partons are formed. The production rate of such events with a gap 
of at least two units of rapidity is estimated [|l(J to be about 10~ 5 of all hadronic Z° 
decays, but with an unknown background from fluctuations in Z° — ► qq fragmentation. 

Here we present the first measurements of rapidity gaps in e + e~ — ► hadrons, using 
Z° decays. We compare our results with the JETSET 7.4 [IT| jet fragmentation model 
and with the perturbative QCD prediction jnj for the rate of large-gap events. We 
compare our results with measurements from pp collisions, where hadronic jets are 
initiated by scattering of u and d quarks or gluons, and with measurements from ep 
collisions, where jets are initiated predominantly by u or d quarks. We have isolated 
event samples enriched in primary light (u,d,s) and heavy (b) quarks, and have mea- 
sured rapidity gap spectra in both. Finally, we have studied the dependence of the 
rapidity gap spectra on the event jet topology. 

The e + e~ annihilation events produced at the Z° resonance by the SLAC Linear 
Collider (SLC) have been recorded using the SLC Large Detector (SLD) [|T2J. This 
analysis used the charged tracks measured in the central drift chamber |TJ[ and in 
the vertex detector WM. The trigger and initial hadronic event selection criteria are 
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described elsewhere ||15| . Events were required to have a minimum of 5 well- measured 
tracks [15|j, a thrust axis [[16] direction within |cos#t| < 0.71, and a charged energy 
calculated from the selected tracks assuming the mass of at least 20 GeV. From 
our 1993-95 data samples 101,676 events passed these cuts, including a background of 
0.3 ± 0.1% dominated by Z° — > t + t~ events. 

Particle rapidity y = 0.5 hn(E + p\\)/(E —p\\), where E is the particle energy and p\\ 
its momentum component along the thrust axis of the event, was calculated from the 
measured momentum assuming the 7T ± mass. The charged track rapidity spectrum is 
shown in Fig. 1(a). We ordered the N tracks in each event by their rapidity, which 
defined N — 1 rapidity gaps, Ay, between pairs of adjacent tracks. For each event 
we considered the largest gap, Ay max , the average gap, < Ay > , and the size of 
the central gap, Ay cent , i.e. that gap containing the mean rapidity of the tracks. The 
Ay max distribution is shown in Fig. 1(b); as Ay max increases it rises to a peak around 
Ay max = 1, and for 1 < Ay max < 6 it falls approximately exponentially; for Ay max > 
6 a 'shoulder' is apparent, suggesting that the data sample contains a class of events 
characterized by large rapidity gaps. 

This feature was found to be due to the 0.3% contamination of r + r~ events in the 



qq sample. The y and Ay max distributions for a sample of JETSET 7.4 |]TTJ qq and 
KORALZ 3.8 [I7| t + t - Monte Carlo events, combined according to their Standard 
Model fractions of Z° decays, and subjected to a simulation of the detector and to 
the same selection cuts as the data, are also shown in Figs. 1(a) and 1(b) respectively; 
the simulation describes the data well. The requirement on the minimum number of 
charged tracks per event n c h was increased from 5 to 7, yielding 100,964 events with 
a reduced contribution from r + r~ events of 0.13%. The resulting Ay max distributions 
for the data and simulated events are shown in Fig. 1(b); for Ay max > 1.0 the data fall 
exponentially and are well modelled by the simulation. Similar results were obtained for 
< Ay > and Ay cent (not shown). All further analysis was based upon the requirement 
n ch > 7. 
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Also shown in Fig. 1(b) is the Ay max distribution from JETSET qq events at the 
generator level; all particles, except neutrinos, with lifetimes larger than 3 x 10~ 10 
seconds, as well as 7r°s, were considered stable and were used in the analysis. Compar- 
ison of this with the measured distribution indicates the effects of the detector and of 
the track and event selection criteria, in particular the exclusion of neutral particles, 
on the Ay max spectrum. The generator-level spectrum exhibits the same structure as 
the data, and it is shifted to smaller values of Ay max and has a steeper exponential 
decrease due to the larger number of particles considered in each event. 

We quantified our data sample in terms of the fraction of events f(y a ) containing no 
charged particles in the interval \y\ < yo/2, shown in Fig. 2; it falls exponentially. The 
JETSET qq simulation falls below the data for y a > 3, but the data are well described 
by the simulation including t + t~ events, demonstrating that even for events with n c h 
> 7 the influence of t + t~ event contamination is noticeable at large gap values. In 



Fig. 2 we compare f(yo) with the perturbative QCD prediction |T0[ for the fraction of 
Z° decays into qqqq and qqgg, where in each case two color-singlet parton-pairs are 
produced with rapidity separation y Q . The prediction, for a s (M§) = 0.120 ± 0.008 
|15| |, has a similar exponential falloff as the data, but is between two and three orders 
of magnitude smaller. The effect of neglect of neutral particles (Fig. 1) is much smaller 
than this difference. These results imply that the production of events with large 
rapidity gaps in hadronic Z° decays is dominated by fluctuations in color-exchange 
jet fragmentation. 

We also compare f(yo) with similar measurements from pp and ep collisions, al- 
though, due to differences in selected event topology, particle selection, and definition 
of rapidity, such a comparison can be made in qualitative terms only. The DO Collabo- 
ration has studied || the fraction of events that have no tagged particles in pseudora- 
pidity space between the two highest-^r jets in their multijet event data sample. Jets 
were defined using a cone algorithm and required to have E T > 30 GeV; y Q was 
taken to be the difference in pseudorapidity between the edges of the two jet cones, and 
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a tagged particle was defined to be a tower in the DO electromagnetic calorimeter with 
Et > 200 MeV. For y a < 1 the DO results decrease exponentially (Fig. 2); for y Q > 1 
they indicate a plateau, with f(y > 3) = 0.0053 ± 0.0009 [|. The CDF Collaboration 
has reported similar results, /(0.8 < yo < 4) = 0.00851q;oqi7 0, based upon absence of 
charged tracks of momentum larger than 400 MeV/c between jets. The ZEUS Collabo- 
ration has studied |7j photoproduced events containing two jets with Et > 6 GeV. The 
fraction of events that have no calorimeter clusters with Et > 250 MeV between the 
two jets, defined using a similar cone algorithm as DO but with a larger cone radius, 
decreases with y (Fig. 2) but may reach a plateau; f(y = 1.7) = O.llj^os 0- 

For y < 1 our measurement of an exponentially falling rate of gap events is qual- 
itatively similar to both the pp and ep cases at low y . For y > 1 our measured 
exponential decrease, which is well described by the JETSET model of color exchange 
between the primary quark and antiquark, contrasts with the onset of a plateau in 
the pp and ep data samples, and supports the interpretation of the large-gap events in 
pp and ep collisions in terms of hard colorless exchange processes || £|, |J. 



The JADE jet-finding algorithm [D| was used to define jets and the rapidity and 
gap distributions were studied as a function of the jet multiplicity of the events. We 
considered 2-jet and >3-jet events defined for values of the scaled jet-pair invariant 
mass in the range 0.005 < y cu t < 0.13. We show results for the 12,394 2-jet events at 
Vcut = 0.005 and the 6,885 >3-jet events at y cut = 0.13. The former are events with 
two narrow back-to-back jets, while the latter events contain at least 3 well-separated 
jets where the additional jet(s) are typically due to hard gluon radiation. 

The 2-jet rapidity spectrum (Fig. 3(a)) peaks at \y\ ~ 2.5 and the region \y\ < 2.5 
is depleted but non-zero, reflecting the color field between the q and q. The >3-jet 
rapidity spectrum peaks at a lower value, which is expected from the lower jet energies, 
the jet topology, and higher event track multiplicity. Corresponding effects are visible 
in the Ay max distributions (Fig. 3(b)); the 2-jet events peak at a gap of 1.5 units, and 
there is a long tail at high gap values; the 3-jet events peak at a gap size of about 
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0.7 units and there are few events with Ay max > 2.5. Similar results were observed 
at other values of y cu t and for the < Ay > and Ay cent distributions (not shown). In 
all cases the high-gap tails demonstrate exponential falloff and are well described by 
hadronic Monte Carlo samples subjected to the same selection criteria. 

The rapidity and gap spectra were studied as a function of event primary flavor 
using impact parameters of charged tracks measured in the vertex detector to tag 
samples enriched in light (Z° — > uu, dd or ss) and heavy (Z° — ► bb) flavor. The 65,243 
events containing no track with normalized transverse impact parameter with respect 
to the interaction point b/a^ > 3 were assigned to the light-flavor sample. The 13,269 
events containing three or more tracks with bjo^ > 3 were assigned to the heavy-flavor 
sample. The light-flavor content of the light sample was estimated to be 85% and the b 
flavor content of the heavy sample was estimated to be 89%. Details of flavor tagging 
are discussed in Ref. [fL9] . 

The rapidity spectrum of the b-tagged sample (Fig. 3(c)) is relatively flat out to a 
\y\ ~ 2.4 and falls sharply thereafter. The light quark spectrum peaks at \y\ ~ 0.6 and 
falls more slowly at high rapidity. This difference can be explained by the fact that 
the rapidity of B-hadrons in Z° decays peaks strongly at \y\ « 2.3, corresponding to 
the average B momentum of 32 GeV/c ||20]| , and the large number of B-hadron decay 
products contributes only a small additional spread about this value; additional tracks 
from fragmentation following the B-hadron formation are restricted to lower momenta 
and rapidities. For Ay max (Fig. 3(d)), as well as < Ay > and Ay cent (not shown), 
the distributions for the two flavor-tagged samples are found to peak at roughly the 
same value, but the b-tagged distribution falls faster at high values. These features are 
expected since the rapidity of B-hadrons in Z° decays is limited to 2.8, and their high 
decay multiplicity populates the region \y\ < 2.8 independent of the fragmentation 
process. The behavior of f(yo) is similar for both flavor-tagged samples (not shown), 
although the slope of the exponential decrease is greater in the b-tagged case. 

In conclusion, we have studied distributions of rapidity and of rapidity gaps using 
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charged particles in e + e~ Z° — > qq decays, and their dependence on the jet topology 
and quark flavor of the events. After accounting for backgrounds from t + t~ events, 
we find that the inclusive distribution of the maximum rapidity gap per event peaks at 
about 1 unit and subsequently falls exponentially with increasing gap size. In contrast 
to results from pp and ep collisions we find no evidence for any anomalous class of 
events characterized by large rapidity gaps, and thus support the interpretation of the 
large-gap pp and ep events in terms of exchange of color-singlet objects. Differences in 
peak position and exponent of the falloff between samples of different jet multiplicity 
or flavor can be explained by differences in event topology, charged multiplicity, and 
the effects of B-hadron production and decay. 

We thank the personnel of the SLAC accelerator department and the technical 
staffs of our collaborating institutions for the successful operation of the SLC and the 
SLD. We thank J.D. Bjorken, S.J. Brodsky and B. May for helpful discussions. 



11 



References 



[1] R.D. Field, R.P. Feynman, Nucl. Phys. B136 (1978) 1. 

[2] J. D. Bjorken, Phys. Rev. D47 (1993) 101. 

[3] DO Collab., S, Abachi et al, Phys. Rev. Lett. 72 (1994) 2332. 

[4] CDF Collab., F. Abe et al, Phys. Rev. Lett. 74 (1995) 855. 

[5] DO Collab., S, Abachi et al, Phys. Rev. Lett. 76 (1996) 734. 

[6] HI Collab., T. Ahmed et al, Nucl. Phys. B429 (1994) 477. 
ZEUS Collab., M. Derrick et al, Phys. Lett. B332 (1994) 228. 

[7] ZEUS Collab., M. Derrick et al, Phys. Lett. B369 (1996) 55. 

[8] See eg: TASSO Collab., M. Althoff et al, Z. Phys. C22 (1984) 307. 

[9] TASSO Collab., W. Braunschweig et al, Phys. Lett. B231 (1989) 548. 
DELPHI Collab., P. Abreu et al, Phys. Lett. B247 (1990) 137. 
OPAL Collab., M.Z. Akrawy et al, Phys. Lett. B262 (1991) 351. 
ALEPH Collab., D. Decamp et al, Z. Phys. C53 (1992) 21. 

[10] J. D. Bjorken, S.J. Brodsky, H.-J. Lu, Phys. Lett. B286 (1992) 153. 

[11] T. Sj6strand, CERN-TH.7112/93 (1993). 

[12] SLD Design Report, SLAC Report 273 (1984). 

[13] SLD Collab., M.D. Hildreth et al, IEEE Trans. Nucl. Sci. 42 (1994) 451. 

[14] C. J. S. Damerell et al, Nucl. Inst. Meth. A288 (1990) 288. 

[15] SLD Collab., K. Abe et al, Phys. Rev. D51 (1995) 962. 

[16] E. Farhi, Phys. Rev. Lett. 39 (1977) 1587. 

[17] S. Jadach, B.F.L. Ward, Z. Was, Comp. Phys. Comm. 66 (1991) 276. 

[18] JADE Collab., W. Bartel et al, Z. Phys. C33 (1986) 23. 

[19] SLD Collab., K. Abe et al, SLAC-PUB 6569, to appear in Phys. Rev. D 
1996). 

12 



[20] T. Behnke, Proc. XXVI International Conference on High Energy Physics, 1992, 
ed. J. R. Sanford (AIP 1992) 859. 



Figure captions 

Fig. 1: (a) Normalized distribution of the absolute rapidity of charged particles; data 
(solid circles) and Monte Carlo simulation (histogram), (b) Normalized distribution of 
event maximum rapidity gap. For n c h > 5 (n c h > 7) the data are shown as open (solid) 
circles and the combined qq and t + t~ Monte Carlo sample is shown as the dot-dashed 
(dashed) histogram. The solid histogram is the distribution for Monte Carlo qq events 
at the generator level. 

Fig. 2: Fraction of events containing no charged particle within \y\ < yo/2 (solid cir- 
cles). Prediction from JETSET qq events (dashed histogram), and from the combined 
sample of qq and t + t~ events (solid histogram). The hatched region is the prediction 
from Ref. |I0| . Similar results are shown from DO [[| (squares) and ZEUS (open 



circles) . 

Fig. 3: Normalized distributions of (a) charged particle rapidity and (b) event max- 
imum rapidity gap for 2-jet (solid circles) and >3-jet (open circles) events; the sim- 
ulations are shown as solid (2-jet events) and dashed (> 3-jet events) histograms. 
Normalized distributions of (c) charged particle rapidity and (d) event maximum ra- 
pidity gap for u,d,s- (solid circles) and b-tagged (open circles) samples; the simulations 
are shown as solid (u,d,s-tagged) and dashed (b-tagged) histograms. 
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